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Abstract 

I show that when non-hnearities are taken into account the Landau theory of 
Fermi hquids predicts the existence of hyperbohc waves in fermionic systems. 
The zero sound is described by a infinite set of coupled non-hnear partial 
differential equations, one for each harmonic of oscillation of the Fermi surface. 
When just a few harmonics are included it is possible to describe the dynamics 
of the first sound in a very simple way. These results lead to the interesting 
experimental possibility of having non-linear effects in quantum liquids. In 
particular I discuss the problem of generation of second harmonic in a Fermi 
system. 

PACS numbers: 47.35. +i, 47.20.Ky, 67.55.Fa, 71.27.+a 
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It is already well established that the dynamics of sound waves in neutral systems such 
as He^ or plasmons in charged systems can be well described by the Landau theory of Fermi 
liquids [jl|. The usual approach is based on the linearized Boltzmann equation from which 
response functions and the dispersion of the collective modes are obtained. It is known 
that the same type of linearized equation leads to the correct description of the low energy 
excitations of the Tomonaga-Luttinger model in one dimension P]. This result means that 
in one dimension, due to the lack of a particle-hole continuum at low energies and long 
wavelengths, the bosonic excitations are simply sound waves. It has been shown recently by 
the generalization of the bosonization procedure in higher dimensions that the actual bosonic 
excitations of fermionic systems at low energies are the displacements of the Fermi surface, 
that is, particle- hole continuum plus collective modes 0. Thus although the quasiparticle 
residue is zero in an one dimensional model with interactions, the Landau expansion for the 
energy and the Boltzmann equation are still valid. 

Amazingly enough, the Landau theory is capable of predicting the existence of a first 
order phase transition. This transition happens at the point where the compressibility 
vanishes (spinodal line) and the system becomes unstable to long wavelength fluctuations. 
The Landau theory predicts that there is an exponential growth of the collective mode inside 
the many-body system leading to phase separation. This is also known as Pomeranchuk's 
instability @] and has been used to study phase separation in mixtures of helium, processes 
involving heavy ion collisions at high energies and more recently, the problem of phase 
separation in two dimensional systems such as some cuprates and nicklates 

The main result of the linearization of the Boltzmann equation is that it leads to plane 
waves which propagate in space with a dispersion defined by the poles of the density-density 
response function In the case of instabilities the frequency of oscillation of the wave is 
found to be imaginary leading to an exponential growth of the amplitude 0. In principle, 
for short range interactions, the growth in unbounded. However, in real systems we always 
have long range forces (such as the Coulomb force) or dissipative mechanisms which tend 
to inhibit the growth of the wave leading to saturation of the wave amplitude 0. 
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I am going to show that, within the framework of Landau theory, if non-hnear terms 
are taking into account, the physics of the waves which propagate in Fermi systems is not 
described by ordinary plane waves but by the so-called hyperbolic waves [||. The non- 
linearities appear when we consider fluctuations of the Fermi surface which are comparable 
to the Fermi momentum. This is similar to the problem of the inclusion of large momen- 
tum scattering in one dimension (backscattering or umklapp) which leads to Sine-Gordon 
theories f^. The terms that we are going to keep shall lead to non-linear bosonic theories 
in dimensions higher than one f^. 

I am going to review the main points on the physics of waves in Landau theory. In this 
way it will become clear the origin of the non-linearities I am talking about. The starting 
point is the Landau theory which assumes that the local change in energy due to a local 
change in the quasiparticle distribution, 6n{k,r,t), is given by , 



where e° is the bare dispersion relation, U (r , t) is an external potential applied to the system, 
V is the volume and p is the quasiparticle two-body interaction (for simplicity I am going 
to assume that these interactions are homogeneous, that is, do not depend on the position 
of the quasiparticles) . 

Then the actual quasiparticle dispersion is obtained from (1), 



It will become clear that non-linearities appear because the quasiparticle dispersion depends 
on the quasiparticle distribution. This is equivalent to problem of flow of fluids where the 
velocity of propagation depends on the density. At this level the equations of motion for the 
quasiparticles are derived from (2) as. 






(2) 



The main problem now is how to describe the evolution in time of the quasiparticle 
distribution given the equations (3). In the Landau theory the quasiparticle distribution, 
n{k,r,t), changes because the scattering of quasiparticles and it is described by the Boltz- 
mann equation, 

where I{n) is the collision integral. Eq.(4) is written in terms of the partial derivatives with 
help of (3) as, 

^I^'^lA + v,-^(fc, r, t) ■ V^e,^{f, t) - Vf:n{k, f, t) ■ V.-e,^(f, t) = I{n). (5) 

As usual we expand the distribution around its equilibrium value, that is, 

n{k,r,t) = UQ^k) + 6n{k,r,t) (6) 

where nQ{k) = 0(/i — e^) with fi the chemical potential. Substituting (2) and (6) in (5) and 
using that V^e° = Vk and Vj^no{k) = —Vk6{fi — e°) we easily find, 

^^''^^^'''^^ + ■ V,.6n{k, f, t) + 6{fi - el) I ■ V,-.f/(f, t) + yY. hS' ' Vr5^(fc', r, t) 

k' 

- ^ E hk' ^k^n[t r, t) ■ \/,5nik', f, t) = I{n). (7) 

k' 

I have kept all terms in the equation and I am not assuming that the deviations from 
equilibrium are small. Observe that the last two terms on the l.h.s. of (7) are not present 
in the usual approach of the Fermi liquid theory because they are second order in the 



deviations. These terms give rise to non-linear effects as I will show. 

Since it is expected that the main physics of the problem depends only on the dynamics 
at the Fermi surface (which is usually the case in fermionic systems in any number of 
dimensions) the solution of (7) can be written as, 

6n{k,f,t) = 6{fi-el)p{k,f,t) (8) 



where p{k, r, t) obeys the following equation, 

k' 

+ ^ E {^kkrk' ■ ^rP{k, f, t)) p{k', f, t) 

k' 

+ ^ E {"^kkk' ■ ^rp{k\ f, t)) p{k, f, t) = i{p), (9) 
k' 

where I have used Pfe/fc^fc/V^5(/i-e^) = -^{.^-(^1)^ k^Pkfkfi) ^"^^ defined E^, = E^/ <^(Ai-e°')- 
For simplicity I assume a spherical Fermi surface. In this way I can simplify further the 
problem by expanding p^ in Legendre polynomials in three dimensions (a completely similar 
procedure can be done in two dimensions in terms of Fourier series 0), 

oo 

p{e,r,t)=VFY.Mr.t) PiicosO) (10) 

1=0 

where 6 is the angle on the Fermi surface that parameterizes k and vp is the Fermi velocity. 
Analogously the quasiparticle interaction can also be expanded as, 

oo 

%'=E/^^KC0S%^,-,). (11) 
1=0 

Since in this paper I will consider the problem of instabilities and non-linearities in Fermi 
liquid theory I will keep only the two first harmonics in the interaction. Firstly it can be 
seen from (9) that in order to have non-linearities in the problem it is necessary to consider 
an interaction which has dispersion at the Fermi surface (V^/^^, 7^ 0) that is, at least fi 
must be different from zero. Also, due to Galilean invariance lll|,p!0|, the effective mass of the 



quasiparticles depends only on /i, namely, ^ = 1 + ^^^J^^ , where A^(0) = y Y^k^ip " ^t) — 
"^^2^ is the density of states at the Fermi surface. Secondly, the compressibility of a Fermi 
liquid, K, depend on the parameter /o [IH (k = „2(]^^j^(o)jQ) where n is the density) and the 
instabilities to phase separation, or Pomeranchuk's instabilities, appear for A^(0)/o < — 1 
when the systems crosses the spinodal line ||^. Thus, in this approach it will be sufficient to 
consider only the first two harmonics, that is. 
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F{k,k')^Fo + FiCOs{ej;g) (12) 



where Fg^,-, = Ar(0)4,-,. 

Notice that in eq. (9) we just need the following summation, 

^ E fU'Pk' = (^o^io + cos , (13) 

k' 

which will be used later in the paper. 

In momentum space the Fermi velocity and momentum are radial and written as, Vk = 
vpr and k — kpr, where vp and kp are the Fermi velocity and momentum, respectively. 
For simplicity I will assume that the external potential depends only on the z direction in 
real space, that is, U{f,t) = U{z,t). Thus it is natural to choose a reference frame in real 
space where the axis coincide with the axis in the momentum space, that is, f = cos 9z + 
sin 9 cos (px + sin 9 sin (py, 9 — — sin 9z + cos 9 cos (px + cos 9 sin (py and (f) — — sin (px + cos (py 
where is the azimuthal angle. Thus the gradient in momentum space is simply Vj^ = 
■^9-^ + fc^ging ^- With this choice the propagation is one dimensional. Of course, more 
general frames can be chosen but this simple choice will help us to understand the dynamics 
of this system. Using the above expressions, from (9) and (13), one finds, 

dp{9,z,t) f dU{z,t) dp{9,z,t) \ 



duo{z,t) 1 ^ dui{z,t) 

Fo cos9 — h -Fi cos 9 — 

oz 3 oz 



+^ sin^ 9iMz^t)pi9,z,t)) ^ 
Skp oz 

which is the final set of equations which must be solved in order to get the dynamics of the 
system. Observe that the non-linearity (last term in the l.h.s. of (14)) comes with higher 
powers in kp^ which shows that this term is important for fluctuations out of equilibrium. It 
is also easy to see that by projecting the above equation in each of its Legendre components 
the final set of coupled partial differential equations has the form {U = 0), 

^^^^ + E cuM)^^^^^ = i{u) (15) 

m=0 



where Cnmiiu}) is a velocity matrix which depends on the components «„. This set describes 
a system of coupled hyperbolic waves (and, in particular, shock waves) where the damping 
depends now on the collision integral I{u). Notice that this is a general result of the Landau 
theory which is obtained just from the fact that the quasiparticle dispersion depends on the 
quasiparticle distribution, eq.(2). 

In order to understand and simplify the problem I consider the hydrodynamic limit where 
the collisions dominate, that is, when it is possible to have local thermal equilibrium in the 
system. In this case the sound mode that propagates is the first sound The dispersion 
of the sound mode can be obtained in Landau theory straight from the compressibility, k 
|l|J^, and it has the dispersion. 



\ f -vf- (16) 

The first sound propagates like a plane wave and at the spinodal line, that is, when Fq < 
—1, the amplitude of the wave grows exponentially given rise phase separation. In the 
hydrodynamic limit, since the collisions dominate, the quasiparticles and collective modes 
decay. In this case it is sufficient to consider only the two first Fourier components in (10), 
that is, 

p{6,z,t) ^ vp [uo{z,t) +ui{z,t)cos6] . (17) 

It is worth mention that in the collisionless limit, when the sound mode is the zero sound, it 
is necessary to solve the complete set (14), or at least, to calculate the non-/mear response 
function of the system |Tl[] . 

By substitution of (17) in (14) and projecting into the Legendre components we find the 
following set of coupled equations. 



dug ^ 2Fit;jr ^ dup ^ 
dt 9kF ^ dz 



^ 1 + ^ + -TTj — «o 



3 V 3 / 9k 



OZ 



dui 4FiVp dui duo dU ui 



where I have introduced in the place of the colhsion integral the relaxation time r in the 
second equation. This term is absent in the first equation because of the conservation 
of number of particles. Here I am going to assume that this term is independent of the 
position. In some cases we can consider r to be momentum dependent which will lead 
to some interesting physics. Notice that (18) resembles the Navier-Stokes equation for a 
fluid with a frictional term If r~^(r) = uVf where u is the viscosity of the quantum 
fluid then we end up with a viscous term exactly as in the Navier-Stokes case. Already in 
this form, eq.(18) has an incredibly large number of interesting physical applications and 
can be studied in different ways. It would be interesting to apply the methods of critical 
phenomena to understand phase transitions in fermionic systems via renormalization group 



methods [|1^,|1J] and relate these equations to the problem of formation of patterns [0. It 
is also worth to point out that if we use a viscous term in the equations we can define a 
Reynolds number, R = CsL/u where L is some characteristic length of the system and we 
can talk about turbulence in Fermi liquids. 

Observe that for kp oo (18) reduces to to the usual linear problem, 

duo j^vf^/^^ F\\ dui _ ^ 



dt 3 V 3 J dz 

dui duo dU ui 

+ vf{1 + Fo)—- = — (19) 



dt dz dz T 

which can be solved by Fourier transform. 



ui [q, uj) = — — (20) 

— c^q — 

where Cg is the sound velocity, given in (16), and U{q, uj) is the Fourier transform of U{z, t). 
Observe that the poles of the solutions above, that is, the eigenfrequencies of the problem 
are given by. 
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When r ^ oo the spectrum is the one expected for a hydrodynamic sound, eq.(16), there 
is no dissipation and the sound propagates as a plane wave. If r is finite and Fq > — 1 



(k > 0) the the roots are pure imaginary numbers for q > {rvp^ A/3{1 + Fo){l + Fi/3))~^, 
that is, cij^ = —ir^{q) where r^(g) is a positive real number. It means that the plane 
wave disappears exponentially, that is, impurities kill the sound mode which is overdamped. 



For Fq > —1 and q < {rvpy 4/3(1 + Fo)(l + Fi/3))~"'^ the sound mode is underdamped. 
However, if Fq < — 1 (k < 0) the roots are purely imaginary and amplitude of the wave 
grows exponentially in time (in a time scale is (r^(g))~^). That is, impurities alone cannot 
stop the phase separation unless r depends on the momentum q. However with a long range 
Coulomb repulsion one changes Fq ^ Fq + ^""^W"^ and the solution of the wave equation 



grows exponentially only for momentum transfer larger than y "^tfi+Fj^^^^ • That is, there 
is a maximum size for the formation of an unstable droplet [0]. Once the instability develops 
in the system, the set (19) is not valid anymore and it is necessary to use the complete set 
(18) and consider non- linear effects. 

In order to understand these effects I will consider the response of this system to a small 
perturbation U. I am going to write an expansion for uq and Ui as, 

uoiz, t) = uo{z, t) + 5uo{z, t) 

ui{z,t) = ui{z,t) + Sui{z,t) (22) 

where uo{z, t) and ui{z, t) are given in (20) and are the linear response to the apphed external 
potential. Suo{z, t) and 6ui{z, t) are the non-linear response of the system to U. In particular 
we are interested in the case where they are proportional to f/^. By substitution of (22) in 
(18) and using (19) we easily find that, 

d6uo{z,t) _|_ "^F A _|_ dSui{z,t) 2FiVf duo{z,t)ui{z,t) 



dt 3 V 3 / dz 9kF dz 

d6ui{z,t) d6uo{z,t) 2FiVf dul{z,t) 

these equations can be solved by Fourier transform. 
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Suppose, for simplicity, that the external potential is monochromatic with some charac- 
teristic frequency fl and wave vector p, that is, U{k,u) = Uo6{k — p)6{u — Q)). Then we 
have, 



uo{k, u) = ^ — ^ 6{k- p)6{uj - n) 



^^0^^' ^) = Q 2n l^,^ 7 ^ ^(^ - 2P)5(^ - 21^) (24) 



which means that the second harmonic is generated in the system. Thus a second peak 
must appear in an experiment which probes density fluctuations. The strength of this peak 
is proportional to the square of the strength of the imposed external potential. Observe that 
the ratio between the height of the peak of the first harmonic to the height of the second 
harmonic is given by, 

UojP^ ^) ^ m* -m p^VL^ ^ 

Experiments in He^ have dealt only with linear effects which are described in (20) |T^. As 
far as I know, there have been no measurements of these interesting non-linear effects in bulk 
He^. In particular it would be interesting to verify the existence of the second harmonic 
generation in He^. The same effect was observed more than thirty years ago in non- linear 
optical systems |jl7 . 
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